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EVIDENCE FOR A TRIAXIAL MILKY WAY DARK MATTER HALO FROM THE SAGITTARIUS STELLAR 
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ABSTRACT 

Observations of the lengthy tidal streams produced by the destruction of the Sagittarius dwarf 
spheroidal (Sgr dSph) are capable of providing strong constraints on the shape of the Galactic gravi- 
tational potential. However, previous work, based on modeling different stream properties in axisym- 
metric Galactic models has yielded conflicting results: while the angular precession of the Sgr leading 
arm is most consistent with a spherical or slightly oblate halo, the radial velocities of stars in this 
arm are only reproduced by prolate halo models. We demonstrate that this apparent paradox can be 
resolved by instead adopting a triaxial potential. Our new Galactic halo model, which simultaneously 
fits all well-established phase space constraints from the Sgr stream, provides the first conclusive ev- 
idence for, and tentative measurement of, triaxiality in an individual dark matter halo. The Milky 
Way halo within ~ 60 kpc is best characterized by a minor/major axis ratio of the isovelocity con- 
tours c/a w 0.67, intermediate/major axis ratio b/a « 0.83, and triaxiality parameter T ~ 0.56. In 
this model, the minor axis of the dark halo is coincident with the Galactic X axis connecting the 
Sun and the Galactic Center to within ~ 15°, while the major axis also lies in the Galactic plane, 
approximately along the Galactic Y axis. 

Subject headings: Galaxy: halo — Galaxy: structure — Galaxy: kinematics and dynamics 



1. INTRODUCTION 

One of the general predictions of structure formation 
within the prevailing cold dark matter (CDM) paradigm 
is that galaxy-scale dark matter haloes should be de- 
scribed by a triaxial density ellipsoid (e.g., Jing & Suto 
2002; Bailin & Steinmetz 2005; Allgood et al. 2006, 
and references therein). The characteristic axial ra- 
tios of these haloes are typically expected to be far 
from spherical (see, e.g., Bullock et al. 2002; Kuhlen 
et al. 2007; and references therein) with characteristic 
central minor/major axis ratios (in the isovelocity con- 
tours) c/a ~ 0.72 and intermediate/major axis ratios 
b/a ~ 0.78 (Hayashi et al. 2007). Adopting the triaxial- 
ity parameter of Franx et al. (1991), 
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this corresponds to a typical value of T = 0.81. As dis- 
cussed by Hayashi et al. (2007), there can be significant 
variation of the axis ratios with radius and between re- 
alizations of similar mass haloes: individual values for 
c/a range from ~ 0.6 — 0.9 while b/a can take values 
~ 0.6 — 1.0. Despite the near-ubiquitous predictions of 
triaxiality from such simulations, there has hitherto been 
little observational evidence to confirm triaxiality in spe- 
cific individual galaxies: while gravitational lensing (e.g., 
Hoekstra et al. 2004; Mandelbaum et al. 2006; Evans et 
al. 2009) and X-ray observations (e.g., Pointecouteau et 
al. 2005) have indicated that non-sphericity appears to 
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be common in dark matter haloes, such studies are sen- 
sitive only to the integral of the density profile along the 
line of sight and do not provide fully 3-D information. 

The Milky Way Galaxy provides perhaps the best lab- 
oratory for testing predictions of halo sphericity since 
the tidal stream remnants of dwarf satellites orbiting in 
the Galactic halo can be traced in three dimensions and 
provide sensitive probes of the underlying mass distri- 
bution. The tidal tails emanating from the Sagittarius 
dwarf spheroidal galaxy (Sgr dSph) 5 have been used for 
such efforts since shortly after the discovery of the dwarf 
by Ibata et al. (1994). Early efforts to model the MW - 
Sgr system (e.g., Johnston et al. 1995, 1999; Velazquez 
& White 1995; Edelsohn & Elmegreen 1997; Ibata et 
al. 1997; Gomez-Flechoso et al. 1999; Helmi & White 
2001) generally concentrated on fitting the orbit of Sgr 
within an adopted Milky Way potential. In more recent 
years however, compelling wide-field views of the exten- 
sive tidal streams associated with the dwarf have been 
provided by the Two Micron All-Sky Survey (2MASS; 
Majewski et al. 2003) and Sloan Digital Sky Survey 
(SDSS; Belokorov et al. 2006), and permitted much 
more detailed exploration of the underlying shape of the 
Galactic potential to be undertaken by Helmi (2004), 
Martmez-Delgado et al. (2004, 2007), Law, Johnston, 
& Majewski (2005; hereafter LJM05), Fellhauer et al. 
(2006), and Cole et al. (2008). Despite the prediction 
from CDM theory that the Milky Way should be best 
described by a triaxial Galactic potential, the majority 
of models to date (although cf. Gnedin et al. 2005) have 
adopted an axisymmetric framework in which flattening 
is introduced only into the Galactic Z axis (i.e. perpen- 
dicular to the Galactic disk). 

5 We refer the non-expert reader to LJM05 (see particularly their 
Fig. 1) for an extended introduction to the general characteristics 
of the Sgr dwarf system. 
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While these recent simulations agree on many aspects 
of the Sgr orbit, no one model has yet been capable of 
reproducing all of the observational data and these mod- 
els reach different conclusions about the shape of the 
Galactic halo depending on which observational data are 
weighted as primary constraints: Some studies favor a 
mildly oblate halo (e.g., Johnston et al. 2005 [hereafter 
JLM05]; Martmez-Delgado et al. 2007), some an ap- 
proximately spherical halo (e.g., Ibata et al. 2001; Fell- 
hauer et al. 2006), and others a prolate halo (e.g., Helmi 
2004). The crux of this "halo conundrum" (highlighted 
by LJM05, and also discussed by Fellhauer et al. 2006, 
Martmez-Delgado et al. 2007, Newberg et al. 2007, and 
Yanny et al. 2009) is that in an axisymmetric Galactic 
potential it is not possible to simultaneously fit both the 
angular precession and distance/apparent radial veloc- 
ity of stars in the leading Sgr stream as it arcs through 
the North Galactic Cap towards the Galactic anticenter. 
In order to match the angular coordinates of stars in 
the Sgr leading arm JLM05 and Fellhauer et al. (2006) 
required mildly oblate or nearly spherical models (with 
Z-axis flattening q z ~ 0.9 — 1.05 in the contours of the 
gravitational potential) for the Galactic dark halo. Such 
models cause the leading arm to reenter the Galactic disk 
in the vicinity of the Sun (to within ^5 — 10 kpc), in 
contradiction with photometric distance estimates (e.g., 
Newberg et al. 2007) from SDSS. Such a model also pre- 
dicts a significant projection of the orbital velocities of 
leading tidal debris onto the observed line of sight, in con- 
trast to spectroscopic radial velocity data (presented by 
Law et al. 2004, LJM05, and more recently confirmed by 
Yanny et al. 2009). In order to reproduce the radial ve- 
locity and distance trends the halo must be prolate with 
q z ~ 1.25 (as demonstrated by Helmi et al. 2004 and con- 
firmed by LJM05), but this in turn cannot account for 
the observed precession experienced by the tidal debris. 

In this Letter we demonstrate that this "halo conun- 
drum" is a consequence of the near-ubiquitous adoption 
of axisymmetric models for the Galactic potential. By 
using fully triaxial models, similar to those predicted by 
standard CDM theory, it is possible to simultaneously 
satisfy all constraints imposed by the structure of Sgr 
debris on the Milky Way's mass distribution. 

2. MODEL 

We adopt a basic formalism similar to that described 
by LJM05 (and described in greater detail by Law et 
al., in prep.) in which the Milky Way is described by 
a smooth fixed gravitational potential consisting of a 
Miyamoto-Nagai (1975) disk, Hernquist spheroid, and a 
logarithmic halo. The respective contribution of these 
components to the gravitational potential is given by: 

^ <2M disk , . 

$disk = -a -. (2) 

^R 2 + (a + Vz 2 + b 2 ) 2 
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where the various constants C±, C2, C3 are given by 
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This form for the halo potential permits flattening to 
be introduced along the three axes qi, q 2 , q z , where q z 
represents the flattening perpendicular to the Galactic 
Plane, while qi and q 2 are free to rotate in the Galac- 
tic Plane at an angle <j> to a right-handed Galactocen- 
tric X,Y coordinate system 6 . Since it is only the ratios 
between these q that have physical significance, we fix 
172 = 1.0 and explore the effects of varying q\ and q z in 
the range 1.0 - 1.8, with = 0°- 180°. We do not 
consider the range q z ,qi > 1.8 since our formulation of 
the Galactic gravitational potential rapidly becomes un- 
physical for larger values. 

We note that this method introduces the flattening di- 
rectly into the gravitational potential (i.e. the isoveloc- 
ity contours) for ease of calculation and consistency with 
previous studies (e.g., LJM05). We have also considered 
more physical models in which axial flattenings are intro- 
duced into the density profile of an NFW (Navarro et al. 
1996) halo model, and the resulting accelerations com- 
puted using the approximate form of the potential given 
by Lee & Suto (2003). Such models give slightly different 
orbital paths, but the qualitative results discussed below 
remain essentially unchanged. 

We assume that the Sun is located 8.0 kpc from the 
Galactic Center and 28 kpc from the Sgr core (Siegel et 
al. 2007). We fix various constants in Equations [2] - [4] 
based on previous work by LJM05, adopting Mdisk = 
1.0 x 10 n M Q , M sphcrc = 3.4 x 1O 1O M , a = 6.5 kpc, 
b = 0.26 kpc, c = 0.7 kpc, rhaio = 12 kpc. Similarly, the 
position, radial velocity, and instantaneous orbital plane 
of the Sgr dwarf are fixed as discussed in LJM05 (see also 
Law et al., in prep). 

For each combination of q±, q z , and <f> tested we pro- 
ceed as follows; (1) We fixed Uhalo in Equation Q] so that 
the LSR has an orbital velocity «lsr = 220 km s" 1 . (2) 
While the dynamics of the leading arm of Sgr tidal debris 
are strongly dependent on the choice of Galactic poten- 
tial, the trailing arm experiences very little precession 
and its characteristic angular position (Majewski et al. 
2003) and radial velocities (Majewski et al. 2004) can 
be well reproduced in any realistic model of the Galac- 
tic potential given suitable choice of the speed Vt an of 
Sgr perpendicular to the line of sight (as demonstrated 
by LJM05). We therefore adjusted t> tan until the best 
match to the trailing debris data was obtained. (3) Mass- 
less test-particles were integrated along the orbit, and the 
quality of fit x to the observational data calculated. 

The accuracy of such test-particle orbits will neces- 
sarily be limited since actual tidal debris (with a range 
of orbital energies and angular momenta) from massive 
satellites will deviate slightly from these orbits; lead- 
ing/trailing arms will fall inside/outside the orbital path 
respectively (see, e.g., discussion by Johnston et al. 1995, 
1999; Choi et al. 2007; Eyre & Binney 2009). However, 
previous work (e.g., LJM05) has demonstrated that for 

6 I.e., when (f> = 0° , q± is aligned with the Galactic X axis and 
Eqn. [3 reduces to * hak) = v£ alo ln([a;/gi] 2 + \ylqif + [z/q z ] 2 + 
r 2 ) 

' halo'- 
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Sgr the test-particle orbits provide a good indication of 
the general trend of tidal debris (particularly for the ob- 
servational constraints upon which we focus) and suf- 
fice for the first-order estimate of the underlying gravi- 
tational potential presented here. 

3. RESULTS 

We first illustrate the spirit of our results with the spe- 
cific case of 4> = 90° (i.e. qi = q y ). In Figure [T] we plot 
radial velocity (with respect to the Galactic Standard 
of Rest [GSR]) and declination along the section of orbit 
leading the Sgr dwarf for models in which the Milky Way 
is taken to have various values of q y and q z . While the 
trend of radial velocity is strongly affected by the choice 
of q z (i.e. different colors in upper panels of Fig. [1]), q y 
has a comparatively minor effect on the velocities (differ- 
ence between solid/dotted/dashed lines of a given color) 
in the key range a ~ 250° — 100°. In contrast, both q y 
and q z affect the precession of the leading arm as charac- 
terized by its declination 6 for a given a or A Q coordinate. 
As shown in the lower panels of Figure [TJ increasing q z 
shifts the arm to greater S for fixed q y (i.e. changing 
from black, to orange, to magenta curves for a given line 
type in Fig. [1]), and increasing q y decreases 5 for fixed 
q z (i.e. changing from solid, to dotted, to dashed curves 
for a given line color in Fig. [T]). These trends indicate 
that while both q y and q z govern the angular precession 
of the stream, it is primarily q z that affects the distance 
to the leading arm, altering the projection of the orbital 
velocity onto the line of sight. 

In the special case of <fi = 90°, q z and qi are clearly sep- 
arable parameters which may be constrained in turn. In 
general however, q z and q\ are not separable and must be 
constrained in tandem so that the resulting tidal stream 
matches both the apparent radial velocity and the ob- 
served angular precession. In Figure [5] we summarize the 
three most relevant observational constraints: trailing 
arm velocity data from Majewski et al. (2004; open cir- 
cles) , leading arm angular coordinates from Belokurov et 
al. (2006; open trianges), and leading arm radial veloci- 
ties from Law et al. (2004; see also LJM05 and Majewski 
et al. in prep.] open boxes). While there is evidence for a 
"bifurcation" in the Sgr leading stream (Belokurov et al. 
2006), we assume that the centroid of the stream is repre- 
sented by the main, highest surface-brightness 'southern' 
branch (see discussion by Yanny et al. 2009; Law et al. 
in prep). As discussed in fj2j it is possible to fit the trail- 
ing arm velocity data in most reasonable models of the 
Galactic potential with suitable choice of the velocity of 
Sgr along its orbit, and all of our models are designed to 
do so. 

We define a x 2 statistic characterizing the quality of 
fit of the Sgr orbit in our various models of the Galactic 
halo to the angular coordinates and velocity data: 

2 1 ( w orbit[«] - V obs .lcad[i]) 2 
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94, n v , trail = 108, and ns = 17, the 



number of observational data points in the 2MASS lead- 
ing/trailing radial velocity samples and the SDSS survey 
fields along the main branch of the Sgr stream respec- 
tively. We adopt an uncertainty a v = I2km s _1 for each 
radial velocity measurement (this represents a combina- 
tion of observational uncertainty and intrinsic stream 
width; see discussion in Majewski et al. 2004), and 
as = 1.9° for each survey field location from Belokurov 
et al. (2006) 7 . The radial velocity w or bit and angular 
location (S or bit of the orbit at each longitudinal position 
(A[i\ or a[i] respectively) is determined via linear inter- 
polation of the orbital path. 

In Figure[3]we show the values of x resulting from vari- 
ous choices for 0, q\, and q z . Values of 4> ~ 90° — 105° (i.e. 
qi nearly aligned with the Galactic Y axis) are strongly 
favored and are the only models in which values of \ < 5 
can be obtained. Given the limitations of the orbit- fitting 
method, it is not possible to conclusively discriminate 
between q\-q z pairs in the triaxial region of parameter 
space q z ~ 1.2 — 1.4 (green/cyan/blue lines respectively 
in Fig. [3]) and qi « 1.2 — 1.8 since deviations of the actual 
tidal debris from the orbital path of the bound satellite 
are expected to have an effect comparable in magnitude 
to the computed values of \ w 3. More detailed fitting 
incorporating finer sampling scales in q Zl qi, and cf> in 
combination with comprehensive N-body models is be- 
yond the scope of the present work and will be presented 
by Law et al. (in prep); at present we note simply that 
the absolute minimum (\ — 2.7) occurs for cj> = 90°, 
q z = 1.25, qi = q y = 1.5. The Laplacian of this potential 
is everywhere positive, indicating that it corresponds to 
a physically realizable density distribution. 

4. DISCUSSION 

In LJM05 we conducted an extensive parameter space 
search exploring the effects of varying the Galactic 
disk/bulge mass, the dark halo scale length, Sgr kinemat- 
ics, the distance to the Galactic Center, the distance to 
Sgr, and the total mass scale of the Milky Way No vari- 
ation of these parameters gave rise to a model in which 
it was possible to simultaneously reproduce all three ob- 
servational constraints (trailing velocities, leading pre- 
cession, leading velocities) on the Sgr stream in an ax- 
isymmetric Galactic halo. In Figure [2] we overplot on the 
observational constraints the orbits in an axisymmetric 
halo that best fit individual constraints: the black curve 
is chosen to best-fit the angular precession of the lead- 
ing arm (qi = 1.00, q z = 0.97), while the green curve is 
chosen to best fit the radial velocity trend of the leading 
arm (qi = 1.00, q z — 1.25). These models have quality- 
of-fit parameters \ — 8.3/9.4 respectively. In contrast, 
the red curve in Figure [2] represents our best-fit orbit in 
a triaxial halo with q z — 1.25, qi — 1.50, <fi — 90°: in this 
model the orbit of the Sgr dwarf simultaneously repro- 
duces both the observed run of angular precession and 
radial velocity for the leading arm. The corresponding 
ratios for the minor/major (c/a = 0.67) and intermedi- 
ate/major (b/a — 0.83) axis ratios in this triaxial halo 
are broadly consistent with the typical values predicted 
by cosmological simulations (e.g., Hayashi et al. 2007), 

7 This uncertainty is a rough estimate of the uncertainty in the 
centroid of the stream at a given longitude based on the apparent 
width of the stream. 
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Fig. 1. — Model orbits are plotted 1.2 Gyr in advance of the current position of the Sgr dwarf (denoted by filled circles) illustrat- 
ing the trend of declination (<5; bottom panels) and apparent radial velocity (i>gsr; top panels) expected for leading arm tidal debris. 
Black/orangc/magenta curves respectively denote haloes in which q z = 1.0/1.3/1.6 while solid/dotted/dashed curves correspond to values 
q y = 1.0/1.3/1.6 (i.e. <j> = 90°) respectively. Given the split in the literature between describing position along the orbit in terms of the 
right ascension (a) or the more natural orbital longitude (A0; see Majewski et al. 2003) we show these relations against both a and Aq 
in left/right- hand panels respectively to aid comparison to plots published by previous authors. 



although the best-fit triaxiality parameter T 0.56 is 
somewhat lower than anticipated. 

While the fit may change somewhat with detailed N- 
body models, we can conclude not only that the Galac- 
tic dark halo may be triaxial within the orbit of Sgr 
(r < 60 kpc), but that the short/long axes of the halo 
may be aligned with the Galactic X/Y axes respectively 
to within ~ f5°. If the axial flattenings are introduced 
into an NFW model of the dark matter density distri- 
bution, qualitatively similar results are obtained with 
Qp,x/qp,y/q P ,z — 1.0/2.0/1.66. As expected (e.g., Kuhlen 
et al. 2007) the isovelocity contours are more spherical 
than the isodensity contours, but both indicate a similar 
preference for a triaxial system in which the major/minor 
axes are approximately aligned with the Galactic Y/X 
axes respectively. 

This alignment is somewhat similar to that of the stel- 



lar halo described by Newberg & Yanny (2006), who 
found a major axis similarly aligned with the Galactic Y 
axis to within ~ 20 — 40°, albeit with a minor axis along 
Z rather than X. There is no clear correlation however 
with the triaxiality of the Galactic bulge: the long axis 
of the bar is currently though to lie within ~ 15 — 20° of 
the Galactic X axis (e.g., Nakada et al. 1991; Morris & 
Serabyn 1996; Babusiaux & Gilmore 2005), close to the 
minor axis of the dark halo. 

In closing, we note two complications of the above re- 
sults. First, our analysis has assumed that the disk lies 
in one of the symmetry planes of the halo potential and 
that the axes themselves do not twist with radius. While 
such assumptions should be explored with future studies, 
we do not expect the broad results of our study to change 
— an orientation significantly out of a symmetry plane 
would result in warping of the disk, which would dampen 
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Fig. 2. — Model orbits are plotted for three models: Axisymmetric halo best fitting leading arm precession (black curve), axisymmetric 
halo best fitting leading arm velocities (green curve), and triaxial halo model (red curve). Solid/dashed lines indicate the orbital path of 
Sgr leading/trailing its current position for 1.2 Gyr. Open triangles represent SDSS data from Bclokorov et al. (2006), open squares/circles 
represent leading/trailing arm M-giant radial velocity data from Law et al. (2004) and Majewski et al. (2004) respectively. Note that 
although the solid black curve fits the angular coordinates well it provides a poor match to the radial velocities; the solid green curve 
provides a good match to the radial velocities but a poor match to the angular coordinates; but the red curve matches both velocities and 
angular coordinates well. 



the misalignment on timescales of a few disk orbits (e.g., 
Dubinski & Kuijken 1995; Bailin et al. 2007; Jeon et al. 
2009). Second, in the family of models presented here 
the Galactic Z axis is the intermediate axis of the triax- 
ial ellipsoid. In a global sense such models are somewhat 
unsatisfactory since orbits about this axis (i.e. in the 
plane of the Galactic disk) are expected to be unstable 
(e.g., Binney et al. 1981). At present, however, such 
models are the only ones in which it is possible to simul- 
taneously reproduce all of the observed characteristics of 



the Sgr stream. 
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